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2 ■ Abstract. We simulate the control of the spin states in a two-electron double 

quantum dot when an external detuning potential is used for passing the system 
through an anticrossing. The hyperfine coupling of the electron spins with the 
surrounding nuclei causes the anticrossing of the spin states but also the decoherence 
of the spin states. We calculate numerically the singlet-triplet decoherence for different 

C^ ' detuning values and find a good agreement with experimental measurement results of 

a similar setup. We predict an interference effect due to the coupling of Tq and T-|_ 
states. 
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1. Introduction 

The development of a working quantum computer, based on electron spins confined in 
quantum dots [Il|2], has been the objective of a large field of experimentalists and 
theorists. One of the most important steps in realizing a spin qubit is achieving 
rapid control of the spin states. For two-electron quantum dots, a method based on 
electrical control of a singlet-triplet qubit was presented by Petta et al. [3], enabling 
control times of a nanosecond time scale. Recently, the control of a singlet-triplet qubit 



based on dynamic nuclear polarization and inhomogeneous magnetic field has been 
demonstrated ^. 

The electrical control of spin states is based on Landau-Zener transitions between 
singlet and triplet states. These transitions take place when external detuning voltage 
is swept over singlet-triplet crossings. The coupling between singlet and triplet states 
in GaAs is due to the hyperfine interaction between spins of the confined electrons 
and spins of the surrounding nuclei |5l|6]. The transition probability depends on the 
sweeping rate of the detuning voltage and on the strength of the hyperfine coupling 
according to the Landau-Zener formula [714T0]. The external magnetic field splits the 
triplet states into three states separated by the Zeeman splitting, T+, Tq and T_. For 
certain detuning values, S — Tj^^ and S — Tq energy differences are smaller than the 
hyperfine coupling. Hence, transitions from the singlet to triplet states are possible and 
result in a nonvanishing occupation probability of the triplet states. 

Petta et al. pTl[T2] studied experimentally consecutive Landau-Zener transitions 
between S and T+ states. They initialized a two-electron double dot so that the two 
electrons are in a single potential well. The Pauli exclusion principle determines that 
the electrons are in a singlet state. After the initialization, they rapidly moved the 
detuning to a value e, where the potential has a double-well form and the electrons can 
tunnel between the wells. They held the detuning fixed for some time, and then rapidly 
moved the detuning back to its starting value, where they measured the probability to 
return as a singlet. They observed an interference pattern of decaying oscillations of 
the singlet probability as a function of e and time. The decay of the oscillations can be 
explained by the fluctuations of the nuclear spins. This experiment has recently been 
studied theoretically in Ref. [T3] . 

In this paper, we study the electrical control of a two-electron double quantum 
dot using Landau-Zener transitions, following the cycle given above. We calculate 
numerically the dynamics of the setup using a 4 x 4 Hamiltonian matrix, derived by 
Coish and Loss [14J , and evaluate numerically the probabilities of the singlet and triplet 
states as a function of time. We observe a similar interference pattern as observed 
in the experiments [11]. The occupation of the Tq state has to be taken into account 
for smaller detunings. The combination of S" — Tq and 5* — T+ oscillations creates a 
secondary interference pattern, which demands higher precision in order to be observed 
in the measurements. 

2. Model and Method 

We model the two-electron system with the Hamiltonian 

where the effective mass m*=0.067me and permeability e=12.7 are material parameters 
for GaAs. The external potential V is divided into two parts, confinement potential Vc 
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Figure 1. (Color online) (a) The confinement potential Vc with L~40 nm, £=5.3 nm 
(dashed), and £=2.7 nm (solid), (b) The energy difference AE between the singlet state 
iS" and triplet states T-,Tq and T_|_ as a function of the detuning e. In the calculation 
of these energies, the hyperfine field is set to zero. The difference of S and Tq states 
is the exchange energy Jq. T_ and T+ states are separated from the Tq state by the 
Zeeman splitting ez- The dot distance is 80 nm and magnetic field 5=100 mT. S and 
T_|_ states cross at £=4.4 nm. 



and potential due to the Zeeman interaction Vz- The smooth two- minima confinement 
potential, shown in Fig. [TJ^a), is of the form 



Vc'(r) = -m*co'^[x'^ — \tanh{x/5)x — ex] + -m*uj^y'^, 



(2) 



where the confinement strength is hu=3.0 meV and the parameter 6=10 nm causes 
smoothing of the potential around x=0. The detuning voltage is taken into account by 
the parameter e. The A parameter is used to fix the coordinates of the potential minima 
to X = ±L, y = 0, which is obtained using 

A = L6 / [L{1 / cosh^ {L/ 6)) + 6 ta.nh{L/ 6)]. (3) 

We use L=40 nm, and the distance of the dots is 2L=80 nm. The potential caused by 
the Zeeman interaction is 



Vzir.s) =^>ijB(r) • s. 



(4) 



where the Lande factor of GaAs is g*=-OAA. The magnetic field can be divided into 
a homogeneous external magnetic field B^xt and an inhomogeneous random hyperfine 
field B„„c(r)- 

We discretize the Hamiltonian using finite difference method and determine its 
eigenvalues using Lanczos diagonalization [15] . We set a numerical grid on the quantum 
dot, 30 grid points in x-direction and 15 points in y-direction. We simulate the 
fluctuations of the hyperflne fleld, causing the decay of the singlet-triplet oscillations, 
by constructing a random inhomogeneous hyperflne fleld. We evaluate the fluctuations 
of the hyperflne fleld by assigning a random hyperflne fleld vector to each grid point. 

In Fig. Wljo), the energy differences between the eigenenergies of the singlet and 
triplet states are shown as a function of detuning for B= 100 mT. The initial value of 



the detuning is chosen higher than 5.5 nm so that the system is initiahzed in the singlet 
state. When the detuning is lowered, the energies of the triplet states approach the 
energy of the singlet state. Around e= 4.4 nm, S and T+ states become degenerate. 
For lower detunings, the energy of the Tq state approaches asymptotically the singlet 
energy. For detunings lower than e = 5 nm, the singlet-triplet energy differences are 
so small that transitions from singlet to triplet states are possible. The coupling of the 
different spin states is via the hyperfine field term. For e < 3.3 nm, the Tq state is closer 
in energy to S state than T+ state. Thus, the Tq state gains larger probability than T+ 
state in the transition from the S state. We observe that the energy differences are not 
linear as a function of the detuning. A linear approximation would be valid very near 
to the crossing points. Hence, it is necessary to use energy values calculated here. 

The energies of the excited states above the four lowest-lying states are much higher 
and need not to be taken into account in our analysis. We approximate the dynamics 
of the system using a 4 x 4 Hamiltonian, constructed in the basis of the singlet and 
three triplet states. For the integrals needed in the calculation of the elements of the 
Hamiltonian matrix, we use the shorthand notation 

K = 9*fJ'B J J dridr2iJ*s/T{ri, r2)5"(ri)^s/T(ri, ra), (5) 

where the space parts of the wave functions "ips/r^ '4's/t are either the singlet or triplet 
wave functions, depending on the calculated matrix element, i is the index of the confined 
electron (1 or 2), and a is x, y, or z. For brevity, we introduce variables [H] that depend 
on hi^2 

h- = \{K + K). 5h- = \{h1-K). (6) 

h^ = h^ ± ih\ 6h^ = dh"" ± idh^ . (7) 

With these variables, the expression for the Hamiltonian matrix H^jf in the basis 
{S, T+, To, T_} reads [H] 

/ -5h+/^ Sh' Sh-/V2 \ 

-5h~V2 Jo{e) + ez + h' h'/y/2 

5h' h+/y/2 Jo{e) h-/V2 

V 6h+/^/2 h+/V2 Me) - ez - h' J 

where ez=g*fiBB^, and Jo is the energy difference between S and Tq states, shown in 
Fig. Wljo)- As the detuning e depends on time, the effective Hamiltonian is also time- 
dependent. 

In order to numerically model the experiment of Petta et al. [11] and compare the 
obtained singlet probabilities with the experimental data, we write the wave function 
of the system in the basis of the four singlet and triplet states, and the singlet and 
triplet probabilities are calculated from the wave function [16]. The dynamics of the 
wave function is calculated using the relation ip{t + At) = exp{—iAtHeff{t)/h)ilj{t). 

In our simulation, we describe the voltage by a detuning parameter e which has 
length as its unit. For the initialization at t=0, we set detuning 3 nm above its minimum 



H. 



eff\ 



a)i 



1.Ub 
1' 


— »n.H^ 


1 




—•—experiment 
theory 


0.95 


A 




0.9 




1 


l\ k j^ 


0.85 






illN^ 


0.8 






1^ 


0.75 






V 


% 






V 


) 5 


10 


15 20 2 






t 


(ns) 





b) 



0.8 


\ 




h 


-A 


A 


0.6 


\V- 


/ \ 


/ \ 


0.4 


\ , 


' \ 


/ 


/ \' 


0.2 




V 


V 


/ 


Vy 


/ \ 


J 





D.34 
0.33 


/ 


\/Wvw 


0.696 
0.695 




\j\jm 






19.7 


19.8 


19.7 


19.8 





10 

t(ns) 



15 



20 



Figure 2. (Color online) (a) The singlet probability Ps, as a function of the time 
t. The blue spheres correspond to measurement results of Petta et al. (upper part of 
Fig. 3(d) in Ref. [H]), black curve represents our theoretical fit for detuning e—A.O nm, 
averaged over 1000 realizations. The theoretical data is shifted by 10 ns, because there 
is a delay in the experiments before singlet-triplet oscillations emerge, (b) The singlet 
probabilities, as a function of the time t. The detuning e is 4.3 nm. The gray line is 
singlet probability calculated for a single realization, black line is an average over 1000 
realizations. The insets show the probabilities when detuning is changed back to its 
initial value in 0.2 ns. After passing the Landau-Zener crossing during the detuning 
sweeping back to its initial value at 19.8 ns, the oscillations are rapidly damped. 



value, which guarantees that the system is initiahzed in the singlet state. We model the 
sweeping of the detuning by lowering it to a minimum value in 0.2 ns. This minimum 
value varies between 2.5 and 5.5 nm. 

The parameters used in our simulation are determined by fitting the theoretical 
results to the experimental singlet probability shown in the upper part of Fig. 3(d) of 
Ref. [H]. The fitting is shown in Fig. [2t^a). We set the detuning in our fit to £=4.0 
nm so that the oscillation periods, which depend only on the S-T^ energy difference, 
are the same. We omit the first 10 ns of the experimental data when the singlet 
probability stays constant. After 10 ns, the experimental singlet probability drops to 
0.7 and then oscillates with period around 4 ns. The oscillation decays towards value 
0.84. We set the hyperfine field strength and the sweeping time to such values that the 
theoretical singlet-triplet oscillation approaches 0.84 and the oscillation decays similarly 
with the experimental data. Best fit is obtained with uniformly distributed hyperfine 
field values between -7.5 mT and 7.5 mT. The root mean square of the hyperfine field is 
4.3 mT. In the experiments, the rms value has been around 2 mT [17]. The difference 
between these values might be due to the different strengths of the confinement potential 
used in the experiments and numerical model. The resulting hyperfine coupling of the 
singlet and triplet states Sh^ is tens of neV, which is in the same magnitude than the 
experimental value 60 neV. We average the spin dynamics over 1000 random hyperfine 
field realizations. As the sum of uniformly distributed random numbers is Gaussian, 
our simulation corresponds to a Gaussian distributed hyperfine field. We sweep the 
detuning to the minimum value in 0.2 ns. With these parameter values, the decaying 
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oscillations are similar with the experimental data. Although the measured value may 
deviate from real singlet probability, this fit gives the relevant scale for the parameters. 

The decay of the oscillations is due to the finite temperature of the nuclear spins. 
Figure Wljo) shows singlet oscillations calculated for a single realization and an average 
over 1000 realizations. We observe for the single realization a non-decaying oscillation. 
For different realizations, the oscillations have different periods and amplitudes. As 
a result, the averaged oscillations decay and approach a constant value, as the solid 
line in Fig. EJ^b) indicates. The detuning is kept in the minimum value for 0-28 ns. 
The sweeping of the detuning back to the initial value in 0.2 ns does not significantly 
change the singlet probability. If we would model the experiment rigorously, we should 
calculate separately each detuning cycle for different waiting times, sweep the detuning 
back to the original value, and save the singlet probability at the end. However, the 
singlet probability does not significantly change during the return sweep. In the insets of 
Fig. Et^b), the singlet-triplet oscillations during the sweeping back to initial detuning are 
shown for a single realization and for an average over 1000 realizations. The sweeping 
takes place between 19.7 ns and 19.9 ns. When Landau-Zener crossing point is passed 
at 19.8 ns, the singlet-triplet oscillation stops and after that the faster oscillations decay 
very rapidly. In practice, the oscillation freezes at the value it had at the crossing point. 
During the 0.1 ns sweep from the minimum value of the detuning to the crossing point, 
the singlet probability changes less than 0.01. Hence, we may end the simulation just 
before the detuning sweep and approximate the final value of the singlet probability by 
using the value the probability has just at the beginning of the detuning sweep. This 
approximation enables calculation of the final singlet probabilities for each waiting time 
in a single run. We only have to make new runs for different detuning values. We set 
the magnetic field to 100 mT, as in the experiment. 

The analysis of the experiment of Petta et al. could be done using a three-level 
model, as it turns out that the T_ state remains unoccupied. It might even be possible 
to obtain an exact analytical solution for the singlet probability. A similar model was 
used to study a recent experiment [H], where a Josephson qubit was coupled to two two- 
level systems. In this case, the two-level systems are not coupled to each other, making 
the calculation of the time dependence easier. We have the additional difficulty of T4.- 
Tq coupling, and the resulting singlet probability should be averaged over Gaussian 
distribution. As our numerical method is reasonably fast to calculate, an analytical 
approach would not be practical for our study. 

3. Results and Discussion 

Now after we have obtained the necessary parameters, we can vary the detuning e and 
also compare results with the experiment of Petta et al. [H]. The calculated singlet 
probabilities for detuning values between 2.5 and 5.5 nm are shown in Fig. [3] (single 
realization) and Fig. H] (average over 1000 realizations). The single realization case is 
shown because the interference effects are then easier to observe. 
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Figure 3. (Color online) Probability Ps of the singlet state S' as a function of time 
t and detuning e calculated for a single hyperfine field realization. The darker colors 
indicate increasing Pg. The white dashed line indicates one of the interference lines 
line due to the S-T^ oscillation and black spheres show the point in the interference 
lines where the singlet probability is diminished due to the S-Tq oscillation. 
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Figure 4. (Color online) Same as Fig. |3l but here the singlet probability Ps is 
calculated using an average of 1000 hyperfine field realizations and the color scale is 
different. 



We have plotted the singlet probability as a function of time (x-axis) and detuning 
(y-axis). The darker colors indicate increasing singlet probability. For detunings above 
the S-T+ crossing point e=4.4 nm, the singlet probability does not change considerably 
from value 1, making the upper part of the figure dark. The period of the singlet 
probability oscillations depends inversely on the energy difference between S and T 
states. In Fig. [3l we observe that the oscillation period diminishes, when detuning values 
are further from the Landau-Zener crossing and the S-T energy difference increases. At 
the crossing, the amplitude and period of the oscillation are largest. The minima of 
the oscillations form observable red lines on the dark background. One of them is 
accentuated with a white dashed line. The same phenomena are present in Fig. HI but 
the averaging makes the figure more blurred. When the detuning is near the crossing 
point, the singlet oscillation drops from 1 to 0.5 and the later oscillations are rapidly 
damped. This causes the yellow band in the middle of the figure. The oscillations for 
detunings above the crossing point are barely visible. For detunings under the crossing 
point, the oscillations are damped and the bands formed by the minima of the probability 
are not so clear. 

The occupation probability of the Tq state affects the singlet-triplet oscillations 
by damping the oscillations for certain detuning values. The red interference lines are 
then slightly narrower. These narrower places in the lines form a secondary interference 
pattern, one of those is indicated with black spheres in the figure. This pattern is better 
seen in Fig. |3l where the widening and narrowing of the lines is clearly observable. For 
example, for detuning e= 2.9 nm, the interference of the triplet states at 3.5 ns causes 
a dip in the peak of the singlet probability in Fig. |5l This effect is quite small and 
not easy to observe experimentally. In the situation where the energy of the 5* state is 
closer in energy to Tq state than to T+ state, the period is given by the S-Tq energy 
difference. Hence, when these two energy differences are equal at e=3.3 nm, the singlet- 
triplet oscillation period has its smallest value 3.2 ns. The interference lines are not 
symmetric with respect to the detuning of the singlet-triplet crossing point. This is due 
to the asymmetry of the exchange energy around the crossing point as a function of 
detuning. As the exchange energy increases rapidly for detunings above the crossing 
point, the period of the oscillation decreases swiftly and the interference lines disappear. 
The features in the figure above e=4.5 nm would vanish completely if the slope of the 
exchange energy were steeper. 

When we compare our results with experiments of Petta et al. [11], we find that the 
interference patterns of the oscillations are quite similar. The contrast of the interference 
lines is better in the experimental figure. This might be due to the small number of 
detuning sweeps used in the averaging of due to correlations between different sweeps. 
The experimental pattern has an area where the singlet probability remains 1 before the 
oscillation begins. The oscillations start later for higher values of detuning, probably 
due to the smoothed voltage pulse profile used in the experiment. We did not model 
this effect in our calculation. For detunings above the S-T+ crossing, our data shows 
small oscillations that are missing from the experimental data. These are, however, 




Figure 5. (Color online) The singlet and triplet probabilities, as a function of the 
time t. The detuning e is 2.9 nm. The dashed line at 3.5 ns marks the point where 
the occupation of Tq state diminishes the probability of the S state. 



much smaller than in the model of Ref. |13j. Probably the measurement accuracy is 
not sufficiently high for the observation of these small oscillations. At the singlet-triplet 
crossing {e ~ 4.4 nm), the numerical singlet probability is smaller than 0.5 at 5 ns. In 
the experiment this minimum is not seen, but the singlet probability goes to 0.5 without 
having decaying oscillations. This might be due to the smoothed pulse proffie, which 
causes over 15 ns delay of the singlet-triplet oscillation at the S-T+ crossing. 

In the measurements of Petta et al. [11], the interference effects due to the 
occupation of Tq state are not visible. The resolution of the experimental data is not 
high enough to observe the interference clearly, but if we compare the data with Fig. |3l 
we find that the widening and narrowing of the bands could have taken place in the 
experiments. With the experimental resolution, the oscillations for the lower detunings 
look quite similar. If the singlet probability could be measured more accurately, it would 
be possible to detect this interference and observe whether both Tq and T+ states are 
occupied. If it is the case that the experiments do not exhibit this interference, the 
reason for the disappearance of S-Tq oscillations could be that the S-Tq coupling is 
much stronger than S-T^ coupling, and the oscillations decay very rapidly. Or the S-Tq 
coupling is much weaker and the S-Tq oscillation period is longer than the measurement 
time. 

A small nuclear polarization (~0.5%) could prolong the singlet-triplet dephasing 
time by two orders of magnitude, as predicted by Ramon and Hu [19]. Then, the 
polarization would not change the S-T+ crossing point but could weaken the S-Tq 
coupling so that the Tq state does not become occupied during the measurement, even 
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for detunings in the vicinity of S-Tq crossing. In that case, the interference effects due 
to Tq state would vanish and the smallest possible period of singlet-triplet oscillations 
would be given by S-T^^^ oscillations at the Zeeman energy, resulting in a minimum 
period 1.6 ns observed by Petta et al. [TT]. 

The results we obtained do not considerably depend on the exact form of the 
potential. We used the potential given by Eq. ([2]), but any other smooth double- 
well potential would produce similar results for the singlet-triplet oscillations. The 
time dependence of the detuning was taken to be linear here. Other monotonic time 
dependencies would not have great effect on the singlet probabilities, because the 
transitions from singlet to triplet states occur at the crossing points, where the time 
dependence can be approximated by a linear fit. The form of the time dependence does 
not matter outside the crossing points. 

In summary, we have studied singlet-triplet decoherence in a two-electron double 
quantum dot with an external detuning voltage. We calculated the singlet probability 
as a function of time for different detunings. The interference pattern of the singlet 
probability is very similar with the recent experimental results and the decay of the 
oscillations is described somewhat better than in another theoretical study of the same 
experiment [13] . Our results indicate that also Tq state could be occupied, even if S-Tq 
crossing point is not within the range of detuning values used. When the system is 
in the proximity of the crossing point, transition to Tq state is possible. In addition, 
we predict an observable effect in the interference pattern due to the occupation of Tq 
state. If the measurement process affects the S-Tq coupling but not the S-T+ coupling, 
the interference effects could disappear, which we suggest might be the case in the 
measurement of Petta et al. [H] . 
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